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We briefly summarize some results obtained by potential NRQCD in the non perturbative regime: 1) a sys- 
tematic procedure to obtain the heavy quarkonium potential in the 1/m expansion (with explicit expressions up 
to 0(l/m ) in terms of Wilson loops), and 2) a description of the heavy quarkonium inclusive decays to light 
particles, and of the NRQCD matrix elements (also the octet ones), in terms of the wave-function at the origen 
and some gluonic field-strength correlators. 



Heavy quarkonium (b-b, c-c) systems have been 
traditionally described by potential models in the 
past, being their inverse size assumed to be of 
0{Kqcd) (and that Aqqd ^ w, being m the 
mass of the heavy quark). Potential models are 
characterized by the introduction of a, more or 
less, phenomenological potential in a Schrodinger 
equation. By assuming some functionality in r 
and by fitting the free parameters of the potential, 
a good description of the heavy quarkonium spec- 
trum was obtained. Nevertheless, there were two 
issues: 1) under which circumstances, and how, 
a pure Schrodinger formulation will emerge from 
QCD in the non-perturbative regime and, if so, 
2) how to obtain the potentials from QCD, or, at 
least, how to relate them with objects eventually 
computable in QCD (then any potential model 
should, at least, be consistent with points 1) and 
2)). The use of effective field theories has helped 
to clarify when point 1) is satisfied and how it can 
be derived from QCD (see [1-3]). The resulting 
Lagrangian of the effective theory (that we name 
potential NRQCD (pNRQCD)) reads 

-CpNRQCD S^(^ido - /l(xi,X2,Pi,P2)js', (1) 

where h is the Hamiltonian of the singlet, i.e. of 
the heavy quarkonium. Actually h is only a func- 
tion of r = Xi— X2, pi, P2, which is analytic in the 
two last operators but contains non-analyticities 
in r. pi = — iVxi P2 = ^iVx2- ^ the 



following expansion up to order l/m^: 

Mxi,X2,Pl,P2) = -^ + -^+W°) (2) 
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The spin-independent terms can be written as 

v^r' - {f ,^p^^^o}+^i^L?+e°^(o,(4) 

and 

where Li = r x pi and L2 = r x p2 and 

V,^^'"\r) - V,^°'^\r;m2 <-> mi). (6) 
The spin-dependent part of ^/(^^o) ^j^g type 

V^f'=vir\r)Li.S,. (7) 

Analogously, for the potential we can write 

V^'jf^ = -Vj^fHr)L2-S2, (8) 

where V^'^g^^r) — V^^'^-'(r;m2 ^ mi). 
For the T/^^'^^ potential we define 

T/(M).^(M)+^(M). 



The spin-indcpcndcnt part can be written as 
V^'/^ = -l{p.-P2,Vy\r)} (10) 

-- ^^|^(Li • L2 + L2 • Li) + y,(i'i)(r), 

while the spin-dependent part contains the fol- 
lowing operators: 

+V^l^'\r)S^-S2 + V^lf{r)S,,{r), (11) 
where Si2(f) = 3f • cti f • cr2 — CTi • cr2 and 

Point 2), i.e. how to obtain the relation of the 
potential with objects computable within QCD, 
can be summarized in the following question: 

How to obtain the potential in terms of Wilson 
loops in the 1/m expansion (also usually named 
adiabatic or Bom-Oppenheimer approximation)? 

The first attempts to answer this question 
started more than twenty years ago. The expres- 
sion for the leading spin-independent potential, of 
0(l/m°), corresponds to the static Wilson loop 
and was derived by Wilson and Susskind [4] : 

V^'^r) = hm^^lniWu). (12) 

Expressions for the leading spin-dependent po- 
tentials in the 1/m expansion, of 0(1 /m^), were 
given in Refs. [5]. The procedure followed in 
these works proved to be very difficult to ex- 
tend beyond these leading-order potentials. In 
Ref. [6], a new method to calculate the poten- 
tials was proposed, where new spin-independent 
(some of them momentum-dependent) potentials 
at 0(1/to^) were obtained. In [7], expressions 
for the spin-dependent potentials were obtained 
in terms of eigenstates of the static limit of the 
NRQCD Hamiltonian in the Coulomb gauge. In 
these works, the potentials did not correctly re- 
produce the ultraviolet behaviour expected from 
perturbative QCD (the hard logs ~ logm). This 
was the first signal that a controlled derivation of 
the potentials from QCD was needed. The solu- 
tion to this problem needs of NRQCD [8] where 



the ultraviolet behavior is encoded in the match- 
ing coefficients of the NRQCD operators. It is 
then possible to incorporate them to the poten- 
tials as done in [9]. At that point, the obtained 
set of potentials at 0{l/mP) seemed to be com- 
plete. 

Nevertheless, this view has recently been chal- 
lenged in Refs. [2,3], where a systematic study 
of the potential has been done within an effective 
field theory framework: pNRQCD. The main im- 
provements achieved in Refs. [2,3] with respect 
these previous computations can be summarized 
as follows: 

A) We have developed a general procedure to 
compute the potential by equating green 
functions in NRQCD and pNRQCD [2]. 

B) We have developed the general method and 
formally provide recursive equations to ob- 
tain the potential at any order in 1/m in 
terms of matrix elements and energies of 
the states solution of the static limit [2,3]. 
These expressions can then be rewritten in 
terms of Wilson loops. 

Points A) and B) solve, in alternative ways, 
question 2) and, thus, finally settle this issue, 
opened since more than twenty years ago. 

Once the formalism has been developed, we 
have been able to obtain, for the first time, the 
complete potential (up to field redefinitions) in 
pure gluedynamics up to 0(1/to) in [2] and up 
to 0(l/m^) in [3]. If compared with the previous 
expressions for the potential that one may find in 
the literature, our novel findings can be summa- 
rized as follows: 

C) A new 0{l/m) potential [2]. 

D) Correct expressions for the spin/ velocity- 
independent terms of the 1/m^ potential: 
Vr^'''^ (K(^'°)) and [3]. 

E) We have clarified that the, so far, standard 
expression obtained by Eichten and Fein- 
berg for vj^^^g^ (^lI'sI) ^^'^ later confirmed 
by several groups was incorrect [3]. 

Let us stress that, to date, A) and B) are the 
only available methods in the literature to com- 



pute the potential in terms of Wilson loops within 
a systematic expansion in 1/m. The attempts 
to implement the method of Eichtcn and Fein- 
berg beyond their leading-order results were not 
able to obtain finite expressions [10]. Indeed, in 
a way, the procedure A) can be seen as the gen- 
eralization of the Eichten and Feinberg method. 
In order to obtain this generalization it was cru- 
cial to understand the computation within an ef- 
fective field theory ideology where equalities be- 
tween Green functions were imposed and interpo- 
lating fields with arbitrary normalizations used. 
The method advocated in Ref. [6] does not ap- 
pear to be correct, at least in its current formu- 
lation, since, for instance, it is not able to ob- 
tain the 1/m potential. The computations in 
Ref. [7] essentially provide the correct expres- 
sions for the spin-dependent potentials (once one 
takes the NRQCD matching coefficients to tree 
level and neglects the tree-level annihilation con- 
tribution in the equal mass case). Nevertheless, 
their methodology needs to be generalized (along 
the lines of [2,3]) to take into account the fact 
that one is dealing with operators instead that 
with numbers in these type of computations. 

Let us also notice that the expressions for the 
potential in terms of Wilson loops encode the full 
perturbative series in a^, at leading order in the 
multipole expansion, in the perturbative situa- 
tion, which is a well defined limit. 

Once h has been obtained, we can obtain the 
energies of the bound states. At the order of in- 
terest, we take the energies from the real part of 
the Schrodinger equation 



(Rc/i) (r|n, /,.s, j) = Enju {r\n,l,s,j), 



(13) 



with quantum numbers n, j, I and s as defined 
in Ref. [11]. Prom the imaginary piece of h, we 

can obtain the inclusive decay widths (to light 
hadrons, leptons or photons) by using the relation 



-2 (n, I, s, j|Im/i|r?,, /, s, j). 



(14) 



This is what we have done in Refs. [12,11] on 
which we report in what follows. 

The imaginary parts of h are inherited from the 
four-fermion matching coefficients of NRQCD: 
Imcfzp Imc4^j, .... Therefore, the structure of 
the imaginary part of the Hamiltonian reads 



Im/. ~ 5:imcf=^ ( A, + + • • • 

k \ 
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being dictated by local potentials. The S-wave 
and P-wave inclusive decay width easily follow 

1 S-wave ~ ^ Im C4_ ^ (15) 
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where RnsOs and Rnjis arc the 5-wave and P- 
wave radial component of {r\n,l,s,j). 

Equations. (15-16) illustrate the factorization 
of the soft and ultrasoft scale since the coefficients 
A]^, Bf-K^^jj, Cj~, etc. are independent of the 
bound-state dynamics. This makes possible to 
obtain model independent predictions from the 
above results, which just come out from the struc- 
ture of our effective theory and the fact that the 
imaginary pieces arc encoded in the four-fcrmion 
matching coefficients. In any case we are able 
to go further and to obtain the coefficients Ak, 
BkA.QQjj, Ck, etc. in terms of gluonic correlators 
or the identity (a specific case of gluonic corre- 
lator). Therefore we can relate them with QCD. 
These correlators could be eventually computed 
in the lattice, evaluated with models or fixed with 
the experiment. For illustration, we show the 
non-perturbative parameter which appears in P- 
wave decays [12]: 



1 /"^ 

C,~— dtt^{gE{t)-gE{0)). (17) 

The explicit expressions for the gluonic correla- 
tors can be found in Rcfs. [12,11], where Eqs. 
(15-16) have been obtained with relative preci- 
sion OiAlco/m'^)- 

The traditional description of the inclusive de- 
cays of heavy quarkonium within an effective 
field theory framework is in terms of expecta- 
tion values of the 4-fermion operators of NRQCD 
[8]. At this respect our procedure can give a 
master equation which relates these matrix el- 
ements with wave-functions at the origen and 
some non-pcrturbative constants independent of 
the mass and the bound state dynamics. Given 
a four-fermion operator O, the master equation 
is the following QH) represents a generic heavy- 
quarkonium state at rest, P = 0): 



{H\O{0)\H) 



1 



(P = 0|P = 0) 
j d^R j d3R'(P = 0|R)(njZs|r) 

(0;xiX2| / ^3^0(0 |Q;x;x^) 



x(R'|P = 0){v'\njls) 



(18) 



where r = Xi — X2, r' = x'^ — Xj, R = (xi -|-X2)/2 
and R' = (x'^ + x^)/2 (note that (R'|P = 0) = 1 
and (P = 0|P = 0) = / d^x). The state 
|0;xiX2) represents the heavy quark-antiquark 
ground state of the QCD Hamiltonian in the 1/m 
expansion (see Ref. [11] for details). 

It is especially important that we can describe 
octet matrix elements, which were thought that 
could not be related with a Schrodinger-like for- 
mulation in any way. We have shown that this is 
not so. Nevertheless, we would like to stress that 
even for the singlet matrix elements, which were 
thought to be related with the wave-function at 
the origen of some Schrodinger-like formulation, 
our results are not trivial because no such for- 
mulation existed in the non-perturbative case be- 
fore. We have provided it, formalizing the con- 
nection between NRQCD matrix elements and 
wave-functions computed in a Schrodinger equa- 
tion and showing how this connection can be im- 
proved in a systematic way. 



Our results allow to reduce significantly, by a 
factor of two, the number of non-perturbative 
constants compared with the number of NRQCD 
matrix elements if we assume that all the states 
below the D-D and B-B production threshold are 
in the dynamical situation Kqcd ^ mtP' . In par- 
ticular, one can think of combinations that are 
free of non-perturbative effects (see [11]). 

I would like to thank N. Brambilla, D. Eiras, J. 
Soto and A. Vairo for collaboration on this topic. 
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